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ABSTRACT: This paper studies gravitational collapse of a complex scalar field at the 
threshold for black hole formation, assuming that the collapse is spherically symmetric and 
continuously self-similar. A new solution of the coupled Einstein-scalar field equations is 
derived, after a small amount of numerical work with ordinary differential equations. The 
universal scaling and echoing behavior discovered by Choptuik in spherically symmetrical 
gravitational collapse appear in a somewhat different form. Properties of the endstate of 
the collapse are derived: The collapse leaves behind an irregular outgoing pulse of scalar 
radiation, with exactly fiat spacetime within it. 
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I. INTRODUCTION 



Recently there has been a lot of new interest in gravitational collapse just at the thresh- 
old for formation of black holes, inspired by the striking numerical results of Choptuik [1] 
on spherically symmetric gravitational collapse of scalar field configurations. Further nu- 
merical results for vacuum relativity in axial symmetry by Abrahams and Evans [2] suggest 
that the phenomena discovered by Choptuik not just restricted to spherical symmetry, and 
a number of extensions have been proposed [3,4,5,6,7]. 

Gravitational collapse has two kinds of possible endstates, according to current views. 
The first endstate consists of a black hole, plus some outgoing matter and outgoing grav- 
itational radiation. The second endstate consists of a stationary remnant star, plus some 
outgoing matter and outgoing gravitational radiation, but no black hole. Choice between 
these endstates depends on initial conditions; gravitational collapse ends in a black hole 
only if the gravitational field becomes strong enough, in the sense of the singularity theo- 
rems [8]. 

A third possibility, the naked singularity, is thought not to occur (cosmic censorship 
conjecture [9]) but is not rigorously ruled out. 

The thought experiment employed by Choptuik, in the context of numerical relativity, 
is to "tune" across the critical threshold in the space of initial conditions that separates 
the non-black-hole endstate from the black hole endstate, and to carefully study the crit- 
ical behavior of various quantities at this threshold [1,10]. Among the infinite number 
of parameters that characterize the initial state, he chooses to tune a single convenient 
parameter p that influences the strength of gravity in the initial state. One could imag- 
ine that the threshold behavior would be intricate, and depend strongly on exactly how 
one tunes across the threshold. However, he finds impressive evidence that the threshold 
behavior is, in important respects, universal. 

For gravitational collapse of fcrmionic matter from suitably regular initial conditions, 
the Chandrasekhar mass sets the scale of minimum mass for any black hole, as long as 
the initial conditions are sufficiently regular. Some condition on regularity of the initial 
conditions is clearly necessary; otherwise one could create black holes of arbitrarily small 
mass with focused beams of ultrarelativistic fermions, evading the threshold set by the 
Chandrasekhar limit. Barring irregular initial conditions, the black hole mass Mbjj(p) 
for fermionic configurations therefore behaves in a simple and rather uninteresting way at 
threshold: It is a step function of p. 

For gravitational collapse of bosonic matter from regular initial conditions, one might 
have thought that Mbjj(p) would likewise be a step function of p, with the role of the 
absent Chandrasekhar mass being played by some mass-scale of the initial data. However, 
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this is not the case. In spherically symmetric collapse of a massless real scalar field 
coupled to gravity, Choptuik finds power-law behavior 

Mbh(p)oc(p-p*)", 7^0.37 (1) 

at threshold, and conjectures that 7 might be universal. For a massless scalar field with 
no self-coupling, no bound star can exist, so that the endstate consists of outgoing scalar 
radiation, plus a possible black hole. 

Furthermore, exactly at the threshold p*, he finds a unique field configuration acting 
as an attr actor for all nearby initial conditions. This field configuration — which we will 
call a "choptuon" — has a discrete self-similarity, by virtue of which it exhibits a striking, 
recurrent "echoing" behavior: it repeats itself at ever- decreasing time- and length-scales 

t' = e-^^t (2a) 
r' = e-"^r (2b) 
ds''^ = e-2"^ds2 (2c) 

cP{t',r')=<Pit,r) (2d) 
(n= 1,2,3,...) (2e) 

where (t, r) are spherical coordinates, and A f» 30 is a constant belonging to the choptuon, 
determined numerically [1]. The solution is thus invariant under a discrete family of scale 
transformations. This choptuon is itself regular, and acts as an attractor for a very wide 
class of regular initial data in spherical symmetry. 

Why study the field configuration exactly at threshold? After all, we understand what 
happens just below threshold (outgoing radiation) and just above (tiny black hole plus 
outgoing radiation). Beyond the obvious fascination of the choptuons, there are at least 
two important motivations: First, choptuons amount to a new kind of counterexample to 
some formulations of the cosmic censorship conjecture, a counterexample that cannot be 
blamed on bad choice of matter fields. Technically they are a counterexample because 
regions of arbitrarily strong curvature are visible to observers at future null infinity. One 
can always bypass this new counterexample by reformulating the conjecture — for instance, 
to specify that generically such behavior does not occur — but the more serious point is 
that choptuons threaten to obstruct any proof of the cosmic censorship conjecture using 
the global theory of nonlinear partial difi^erential equations. Therefore it will be necessary 
to confront and understand them. 

Secondly, and even more importantly, choptuons represent in principle a means by 
which effects of extremely strong field gravity, even quantum gravity, can be observable 
in the present universe. An experimenter, tuning across the black hole threshold in a 
succession of gravitational collapses, should be able observe events in which some outgoing 
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radiation appears from regions where the spacetime curvature is as strong as the Planck 
value, and where quantum behavior of general relativity or string theory may be studied. 
Therefore, the quantum corrections to the classical choptuon — or, indeed, the quantum 
gravitational and stringy generalizations of the choptuon — demand study. 

Striking though they are, the numerical computations involving partial differential 
equations (PDE) do not give complete information about the choptuons. For instance, 
output is restricted to the domain covered by the numerical coordinate system adopted, 
which may not cover the whole domain of dependence of the initial data. This makes it 
hard to study the endstate, the spacetime singularity to which the choptuon collapses, 
and the burst of radiation that it emits. Therefore it is valuable to find further choptuons 
that can be studied by analytic techniques, or mostly-analytic techniques including numer- 
ical solution of ordinary differential equations (ODE). In particular, Evans and Coleman 
[11] obtained choptuons for the spherically symmetric collapse of hot gas, which are con- 
tinuously self-similar, a special case of discretely self-similar choptuons. The Evans and 
Coleman configurations show threshold behavior of Mgjj, but not echoing. Continuously 
self-similar solutions are convenient because they are governed by ODE, not PDE. The 
greater numerical accuracy obtainable with ODEs is desirable, not for its own sake, but 
because of the light it may shed on the deeper questions of universality. 

The subject of this paper is the behavior at threshold in gravitational collapse of a 
complex scalar field (j), under spherical symmetry. The use of a complex scalar field, rather 
than a real one, will allow "echoing" to occur in the form of phase oscillations: 4) changes 
phase, but not amplitude, under a scale transformation. We thereby construct and study 
a continuously self-similar, complex choptuon as solutions of the coupled Einstein-scalar 
field equations. In spherical coordinates (t, r) our scalar field takes the form 

0(t,r) = (-tr/(-r/t), (3) 

where a; is a constant to be fixed by a nonlinear eigenvalue problem arising from the 
field equations. This field is continuously self-similar: under a scale transformation by an 
arbitrary parameter A, it transforms solely by a phase factor: 

t' = erh (4a) 

/ = e^^r (4b) 

ds''^ = e~^^ds^ (4c) 

(j)(t', r') = e-'''^(f){t, r) (4d) 

(0 < A < oo) (4e) 

This solution thus exhibits "linear" phase oscillations — a form of the "echoing" — super- 
imposed on an exact continuous scale symmetry. It is thus similar to, but "less nonlinear" 
than the solution of Choptuik [1] . 
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We are able to construct a choptuon as a nonextendible spacetime, singular only at 
the single point {t, r) — (0, 0), that evolves from regular initial data. A burst of outgoing 
scalar radiation, irregular but of finite energy, emerges just prior the retarded time at which 
the singular point forms at (0, 0), and explodes along the future light cone of that point. 
The most surprising aspect of this choptuon is that spacetime appears to be precisely 
flat throughout the interior of the future light cone of the singular point. Thus, in this 
example anyway, a choptuon leaves behind exactly flat space, with no radiation at all after 
the irregular outgoing burst. 

Since a real ^ is a special case of a complex (f), all of Choptuik's numerical results 
immediately apply as special cases to a complex field; in addition, we find exist intrinsically 
complex choptuons which do not reduce to the real case. Which of these choptuons is the 
strongest attractor is an important question not addressed in this paper; we will return to 
it in future papers. Preliminary numerical evidence by Choptuik [12] seems to show that 
the real choptuon is at least not rendered unstable when complex degrees of freedom are 
added. 



II. FIELD EQUATIONS 
We begin with the spacetime metric in the form used by Choptuik, [1] 

where a{t, r) and a(t, r) are functions of time and radius. This is an example of "radial 
gauge" because the area of spheres (given by the coefficient of dO?') defines the radial 
coordinate. Radial gauge breaks down at an apparent horizon, and so another coordinate 
system will need to be used if an apparent horizon appears. The time coordinate is chosen 
so that gravitational collapse on the axis of spherical symmetry first occurs at t = 0, and 
the metric is regular for t < 0. This metric remains invariant in form under transformations 

t = t{t') (6) 

of the time coordinate, and such a transformation can be used to set a{t, 0) = 1 for t < 
on the axis. By regularity (no cone singularity), a(t, 0) = 1 on the axis for t < as well. 

Matter consists of a free, massless, complex scalar field that obeys the wave equation 

□0 = 

while the Einstein equations are 
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These equations amount to 

1 „ /a 



dr {^r'^drcf) = dt (^9^^) (7a) 

—OrOt Ora — ( 7bj 

a a r 

( —dra + -dro] = Stt { X^drcj)* drcj) + -^dt(j)dt(j)*] (7c) 

— a^a = 47r {dr4>*dt4> + dr(j)dt(j)*) (7d) 
or 

These equations admit global C/(l) symmetries for a constant A, 

(j)' = e^^0, -oo < A < oo (8) 

leaving the metric invariant. 



III. CONTINUOUS SELF-SIMILARITY 



We now derive the form of the scaling transformations, Eqs. (4), for the fields. We 
assume spacetime is spherically symmetric, and admits a homothetic Killing vector field 
[13] C obeying 

^^9liu = ^n^u + ^v^ii = '^Qfiy, (9) 

which generates a continuous one parameter family of homothetic motions (self-similarities) 
on spacetime, Eqs. (4). Here, C denotes the Lie derivative. In a coordinate system, Eq. (5), 
^ may always be taken in the form 

ed,i = tdt + rdr, (10) 

possibly after a further coordinate transformation, Eq. (6). Equations (9, 10) are imple- 
mented by Eqs. (4abc). The metric functions are then of the special form 

a{t, r) = a{-r/t), a{t, r) = a{-r/t). (11) 

The minus sign is chosen so that —r/t > where t < to the past of the singularity. 
Dimensional analysis (in classical general relativity, not quantum field theory) suggests 
that (plt, r) should have dimensions (length)^. Therefore, should we assume that is 
invariant under ^, C^cf) = or (j){t,r) = (f){r/t)? No, the global U{1) symmetry must 
be accommodated too. Quite generally, one must allow for spacetime symmetries to get 
mixed up with internal symmetries; for instance, the Maxwell equations for the potential 
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Afj, are conformally invariant only in the sense that a gauge transformation be allowed 
to accompany each conformal transformation. In this case, we must allow some U{1) 
transformation, Eq. (8), to accompany each scale transformation, 

C^4> = ed^,(t> = iu^<t> (12) 
under an infinitesimal scale transformation, or 

r) — exp(— ia;A)0(t, r) 

under the finite scale transformation, Eqs. (4), with oj a constant of the solution. 

This transformation law for can be conveniently implemented by adopting the form 

0(t,r) = (-tr/(^). (13) 

The time coordinate t has now been redefined, Eq. (6), so that the first singularity of the 
collapse is at t = 0; and a new independent variable z has been introduced by 

z = -r/t; (14) 

z is invariant under scale transformations, Eqs. (4). Also, from Eq. (11), 

a = a{z),a = a{z). (15) 

With the scale symmetry implemented by Eqs. (13, 14, 15), the next step is to solve the 
field equations. We can give initial conditions in the hypersurface t = —1 and then evolve 
forward in time. It follows from a theorem of Berger [14] that the Einstein equations are 
compatible with homothetic symmetry, in the sense that scale invariant initial data will 
always evolve to a scale invariant spacetime; therefore the field equations ought to have a 
solution. 

Remarks should be made as to how suitable the assumption of self-similarity is. Self- 
similar spacetimes can never be asymptotically flat (because the ADM or Bondi mass would 
define a length scale, breaking self-similarity); furthermore, there are some reasons to be- 
lieve that self-similar spacetimes can never be spatially compact [15]. Spatial compactness 
is irrelevant here, but shouldn't gravitational collapse be modelled by an asymptotically 
flat spacetime? However, the self- similar solution should be interpreted as a model of the 
gravitationally collapsing region out to some radius R, where it can match smoothly onto 
a non-self-similar, asymptotically flat region. To anticipate, we will flnd that our solutions 
always have a certain horizon called the past similarity horizon; as long as the matching 
radius R is outside the past similarity horizon, the gravitational collapse that we study 
here will remain entirely within the domain of dependence of the self-similar region of the 
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initial data. For instance, for the complex choptuon discussed below, a value R > 5.004 
suffices. 

Transform the metric variables a, a to a new set b, u given by 

b{z) = a{z)/a{z), 
u{z) = a (z) — 1. 

Under the similarity hypothesis, b = b{z) and u = u{z) are functions of z = —r/t alone. 
Under this transformation, the metric, Eq. (5), becomes 



e2^((l + w) -{b'^ - z'^)dT'^ + 2dTdz + dz'^ +z'^dQ'^^ 



(16) 



where r = Int. Following Choptuik's notation, represent 4> in terms of complex functions 
r),n(t, r)) where 



n 



a 



dt^ = -{-ty^-^p{z) 



a 



where the functions q{z),p{z) are: 



q{z) = 



dz 



p(z) ^ ^ [iujf - z^ 



a 



The field equations now take the form 



d 



Q 



dz\P J z 

db bu 

dz z 
du 

dz 
du 

dz 



f u + 2 
u 



1 



(u+l) 47rz(|gp + - 

L z-i 



— = -87T{u + l)bRe{q*p), 



-z b 
b -z 



u 



P+ 
P- 



(17a) 

(17b) 
(17c) 
(17d) 



where 



P± = iuj + u ± 1, 
A = 6^ - z'^. 

The boundary conditions are as follows. At ^ = 0, regularity of solutions on the axis of 
spherical symmetry demands 

6(0) = 1, 
u(0) = 0, 
q{0) = 0. 
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The boundary value P = p{0) of p is a free boundary condition, and will henceforth be 
taken real and nonnegative by a global phase transformation of 4>, without loss of generality. 
Integration of Eqs. (17) from ^ to ^ = +00 is tantamount to solving the initial value 
problem in the spacelike hypersurface t — —1. Evolution in t is then fixed by self-similarity. 

The wave equation Eq. (17a) has a singular point when A vanishes, i.e., when b{z) — z. 
This value of z will be called Z2, and the point will be called a similarity horizon. It 
represents the null hypersurface in spacetime where the homothetic Killing vector, timelike 
near the axis, becomes null. The value of u there, 

K = 1 - u(z2), 

is free, while b{z2) = Z2 is of course determined by Z2. 

The fiat space wave equation provides a simple introduction; see Fig. 1. Here the metric 

is 

and the general spherically symmetric solution to the massless wave equation is 

/(t^r)+9(t + 

r 

where / and g are arbitrary functions. The homothetic killing vector is ^ = tdt + rdr and 
all self-similar solutions, Eq. (3), are linear combinations of 

0±(t, r) = = {-ty^^t^ (19) 

r z 

where 

z = —r/t. 

Here 0"*" is an outgoing wave, is regular at the past similarity horizon z = +1, and is 
irregular at the future similarity horizon z = —1. In contrast, (f)~ is an ingoing wave, 
is irregular at the past similarity horizon z = +1, and is regular at the future similarity 
horizon z = —1. All solutions are regular at the hypersurface t = or z = 00, except at 
the single spacetime point (t, r) = (0, 0). Only the linear combination (p^ — (p~ is regular 
at the origin r = of spherical coordinates (for all t ^ 0). We can transform from the 
coordinates {t, r) to the coordinates (r, z) where r = Int. Then the metric becomes 

= _ ^2^^^2 ^ 2dTdz + dz'^ + z'^d^^^ (20) 

The coordinates (r, z) have a singularity on the spacelike hyperplane t — Yet another 
coordinate system for fiat spacetime, regular at t = 0, is (p, f), where p = ln(r) and 
V = 1/z = —t/r; the metric is 

ds^ = e^P (^-dv^ - 2dvdp + (1 - v^)dp^ + dn'^^ . (21) 



Returning now to the curved spacetime metric Eq. (5), analysis of the wave equation 
near z = Z2 shows that q{z) and p{z) hkewise have a regular solution and an irregular 
solution at the similarity horizon. The regular solution behaves like 

/ ~ const (22a) 
[QiP] ~ const (22b) 

and the irregular solution behaves like 

f ^{z-Z2f''+^^''' (23a) 
{q,p)^{z-Z2f^+^^'^-^ (23b) 

This behavior is identical to the case of flat spacetime except for one aspect: The presence 
of the quantity 

dh 

K= 1 - — 

aZ Z=Z2 

in the exponent. For flat spacetime k = 1. In a self-similar spacetime, n plays a role 
parallel to that of the surface gravity of the event horizon of a black hole in a stationary 
spacetime, and in this paper we will simply call it the "surface gravity" . 

The regular solution again represents outgoing radiation crossing the similarity horizon, 
and the irregular solution represents ingoing radiation propagating along the horizon. Since 
the subject of study is gravitational collapse from regular initial conditions, and since the 
similarity horizon 2; = 2:2 is in the Cauchy development of the initial data — it is to the 
past of the earliest singular point (t,r)=(0,0) — regularity of will be demanded on the 
similarity horizon. This means that a linear combination of q and p must vanish at 2; = 2:2: 

/3+g(22) - P-p{z2) = 0, (24) 

so that, for instance, we can choose ^(2:2) freely as a complex valued boundary condition 
at the similarity horizon, and qiz'}) is then fixed. 

To sum up, the free data for the system will be taken as 

a;, (25a) 

P = p{0) (real), (25b) 

22, (25c) 

K=l-u{z2), (25d) 

^(22) (complex), (25e) 

amounting to six real constants. Solutions to the system Eq. (17) must be found by 
searching in the six-dimensional data space. 
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IV. CONSTRUCTION AND PROPERTIES 
OF THE COMPLEX CHOPTUON 



The numerical methods we used followed Numerical Recipes [16] Chapter 16: We 
handled the system as a two-point boundary value problem with one fixed boundary (at 
z = 0) and one free boundary (at z = Z2). We "shot" with an adaptive step ODE solver 
from each boundary point, to meet at a point zi in the middle. The six free data values 
were adjusted in an outer loop with a Newton's-method solver for nonlinear equations, 
the six nonlinear equations being the matching conditions for the six ODEs (17) at zi. 
Convergence of the Newton's method then identified a solution on the domain < z < Z2- 
Initial work was done with a Runge-Kutta integrator, while a Bulirsch-Stoer integrator 
was subsequently used for higher precision. The solution was then continued to larger z 
without the need for further boundary conditions. 

To enforce numerically the boundary condition on the horizon is difficult, because the 
irregular part of 0, which we want to annul, vanishes anyway if k < 1, and the smaller 
K, the faster it vanishes. Said diflferently, the diflSculty is that we want to start a purely 
regular solution at the horizon and continue it, but the unwanted irregular solution grows 
rapidly as it is continued away from the solution; this is a numerically unstable situation. 
We handled this by computing analytically the second derivative of the regular solution at 
z = Z2, and using a second-order Taylor expansion to start the solution there. Even so we 
encountered numerical problems for k^1/3, as discussed in Appendix A. 

A single new solution was found on a domain < z < Z2, obtained with the values 

uj = 1.9154446 ± 0.0000001, 

P = p(0) = 0.67217263 ± 0.0000004, 

Z2 = 5.0035380 ± 0.0000002, 

u{z2) = 1-K = 0.39707205 ± 0.00000003, 

p{z2) = (0.020305344 ± 0.000000003) + (0.007153158 ± 0.000000002)z, 

This solution will be called the "complex choptuon." Its complex conjugate solution, with 
the sign of oj changed, also of course exists. Figure 2 displays the functions q{z), p{z)), 
b{z), u{z) for the complex choptuon. 

We also found solution at a; = 0, corresponding to a real scalar field 0; this, however, 
is nothing but the flat {k = 0) Robertson- Walker cosmological model with a real scalar 
field as matter, reversed in time. Further details of the numerical results are reported in 
Appendix A. 

This solution was continued smoothly past 2; = 2:2 to large z. The nature of the point 
2; = 00 on the z-axis must now be clarified. Examination shows this to be a singular 
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point of Eqs. (17). However, on physical grounds there should be no spacetime singularity 
here, since z = oo corresponds to the spacelike hypcrsurface t = 0, which should be regular 
except at the axis, since it lies in the Cauchy development of the initial data. Any apparent 
singularity there must fall under strong suspicion as a coordinate singularity caused by a 
bad choice of time coordinate, Eq. (6), near t = 0. Indeed, the system of equations can be 
rendered regular at ^ = oo by the following change of variables: 



dw 


Z'^ 


w)\ 
w)) 


= z^-'"" 1 


v{w) 
u{w) 


_b{z) 

Z 

= u{z). 



w 



at z = oo. 



qlz 
p(z 



The spacetime metric, Eqs. (5, 16), becomes 

^ g2p ^1^-^ ^ _^^2 ^ 2dwdp + (1 - v'^)dp^ + dn'^^ 



(26a) 
(26b) 

(26c) 
(26d) 

(27) 



where p = Inr. Moreover, the scalar field written as Eq. (3) appears irregular at z, but 
can be written in a regular form as 



'F{w), 



(28) 



and {Q{w), P{w)) are derived in a regular way from F{w). In terms of the new independent 
variable w, the equations of motion become 



_d_ 
dw 



Q 
p 

dv 

dw 
du 

dw 
du 

dw 



1 [ V -1 
l-t;2l-l V 



P+ 
(3- 



u — 1, 
u + 1 



47r(|gp + |P|2) 



u 



-87r(w+ l)Re(g*P). 
From the last two equations, u{w) can be expressed as 



uiw) = 47r[|g|2 + |P|2 + 2vRe{Q*P)]. 



(29a) 

(29b) 
(29c) 
(29d) 

(30) 



The point 2; = oo is now marked by v{w) = 0, and the system is clearly regular at this 
point. (The apparent pole at v = in Eq. (29c) is cancelled by a zero in the numerator, 
from Eq. (30).) The singular points v = ±1 are, in contrast, true singular points of the 
system, and are, respectively, the past similarity horizon at 2; = 2:2, and a new future 
similarity horizon. 
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Integration of Eqs. (29) is tantamount to integration of the field equations in a timefike 
hypersurface r = 1. Due to self- similarity, this in turn is tantamount to evolution in a 
timelike direction of a whole spatial hypersurface < r < oo, in the region outside of any 
similarity horizons. 

Numerically we proceed as follows. Having found a solution for < z < Z2, "i^c integrate 
it a little further in z and then transform by Eqs. (26) to the w variables. At this point v is 
a little less than 1, and decreasing. Then we integrate in w (noting that the solution, if it 
approaches w = 0, it always continues smoothly through) until one of two things happen. 

The first alternative is that u ^ oo at some point; we interpret such a point as an 
apparent horizon in spacetime, where radial gauge coordinates (which we use throughout, 
see Eq. (5)) break down. Such an apparent horizon probably means that the solution 
represents, not a choptuon, but a black hole that is growing by self-similar accretion of 
scalar field. We believe that such a solution, if continued further in a different coordinate 
system, would always encounter a spacetime singularity within the black hole, but we will 
not pursue this issue here. 

The second alternative is that f — > — 1 at some finite w. We interpret this as a second, 
future similarity horizon. The complex choptuon behaves in this way. What are the proper 
boundary conditions at the future similarity horizon? Above, at the past similarity horizon, 
we argued that should remain regular because it arose from regular initial conditions, 
and was not influenced by the spacetime singularity to its future. That argument does not 
hold water here, because the future similarity horizon is outside the domain of dependence 
of the initial data, in particular, it is the Cauchy horizon of that domain. Observers on 
the future similarity horizon will see data coming from the singularity, and indeed it is 
a very interesting question to ask what they will see. Therefore no boundary condition 
at all is enforced on at the future similarity horizon; ((/>, Q, P) may be an arbitrary 
combination of regular and irregular solutions. This is good, because we have already used 
up all our boundary conditions at the axis and the past similarity horizon, and we would be 
embarrassed to have to obey further boundary conditions. However, we must still decide 
how to continue the solution across the horizon, because the horizon represents a singular 
point of Eqs. (29), and furthermore is a Cauchy horizon in spacetime. Evolution across 
a Cauchy horizon in classical general relativity is never unique. Quantum considerations, 
which might fix the evolution in some way, are beyond the scope of this paper, though 
later papers will treat it. Therefore we will continue with a conservative assumption about 
how to evolve. 

Drop self-similarity for a moment, and look at the general boundary conditions for a 
wave equation in curved spacetime. A test wave function is allowed to have discontinuity, 
or a discontinuity in some higher derivative, across the characteristic surfaces of the wave 
equation, which are the null hypersurfaces of spacetime. However, if the wave equation 
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is coupled to gravity, a discontinuity in value of <j) will cause an infinite-mass singularity, 
and is therefore forbidden: T^j^ ~ (V(/)) ~ 5() • However, a discontinuity in the first or 
higher derivative of 4> is allowed. In terms of the characteristic initial value problem, the 
initial data for the wave equation on a characteristic hypersurface is the value of cj) itself, 
but does not include any derivatives of cf). 

We now assume that spacetime continues to be self-similar to the future of the future 
similarity horizon. So turn now back to our self-similar equations for 0. The matching 
rules for at the future similarity horizon are as follows. The wave function (/> must 
be continuous across the horizon, to forbid an infinite-mass singularity at the horizon. 
The regular part of (f) represents incoming radiation crossing the horizon from the regular 
region exterior to it. In contrast, the irregular part of consists of outgoing radiation 
originating very close to the spacetime singularity, streaming along the similarity horizon. 
In w coordinates, the wave variables (Q, P) of the irregular part of the solution oscillate an 
infinite number of times while approaching the future similarity horizon, and also die as a 
power law if k > 0. However, if k = 1 the amplitude of the oscillations remains constant. 
The irregular part of the solution is allowed to be discontinuous; that is, the amplitude of 
the irregular part of (Q, P) can be different on the two sides of the horizon. This causes 
at worst an (allowable) jump in the stress-energy. 

In general we expect to find that (Q, P) contains both regular and irregular parts at 
the future similarity horizon. Then the regular part must be continuous, but the irregular 
part is allowed to jump arbitrarily. The solution can still be fixed uniquely, if we are willing 
to postulate that spacetime is smooth along the future time axis t > 0, r = 0. In that case, 
(j) and the spacetime geometry must obey boundary conditions at the future time axis, 
which are essentially the same as the boundary conditions we have already imposed at the 
past time axis. Just counting degrees of freedom in the boundary conditions, we expect 
this boundary condition to fix a unique solution (or perhaps a discrete set of solutions). 

However, at this point a numerical miracle happens at the future similarity horizon. 
We find that for the complex choptuon, is purely irregular there, so that the regular 
piece of cj) vanishes there, which entails that k = 1 there (to an accuracy of 10-6), so that 
the future similarity horizon carries initial data for fiat spacetime with constant (j). In 
particular the mass aspect vanishes on this null hypersurface. 

Figure 3 displays the functions Q{w), P{w), v{w), u{w). Note that, at the future 
similarity horizon marked by v = —1, u = 0, which means the mass aspect vanishes there; 
this is the numerical miracle discussed above. As indicated in Figure 4, the scalar field is 
irregular at the future similarity horizon, and oscillates and infinite number of times. A 
distant observer, measuring outgoing radiation in the scalar field, would see such a signal 
coming from the threshold gravitational collapse. 
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Consider again the evolution of the solution to the future of the future similarity 
horizon. In view of the boundary conditions, we are allowed to choose an arbitrary irregular 
part just to the future of this horizon. The choice of zero irregular part gives flat spacetime 
to the future of this horizon. Any other choice will give a negative mass in this region, 
and will create a negative mass naked singularity along the future time axis t > 0, r = 0. 
Because we are attempting to evolve across a Cauchy horizon, the choice is not obligatory. 
However, the most sensible choice is clearly the one that gives fiat spacetime. We therefore 
conclude that the complex choptuon leaves behind it flat spacetime to the future of the 
singular point at (t, r) — (0, 0), to an accuracy of 1 part in 10^. 



V. DISCUSSION AND OUTLOOK 



In Figure 4, we show the interpretation arrived at for the complex choptuon. In region 
I, there is a collapsing ball of scalar field, bound by its own self-gravity; it acts as a near 
zone for scalar radiation. The ball collapses toward the spacetime singularity at the origin, 
shown as a single point. As the ball collapses, scalar field is partially trapped by spacetime 
curvature, but also continually leaks out of the gravitationally bound region, across the 
past similarity horizon at v = +1, and radiates outwards through regions Ila and lib. The 
boundary between Region Ila and Region lib is the surface t = in Choptuik coordinates, 
or 2; = oo. However, this is merely a coordinate singularity, and there is no physical 
boundary between the two regions; they should be considered as a single Region II, which 
acts as a transition zone between the near and far field for scalar radiation. Spacetime 
is curved in Region II, and in general some backscatter of outgoing radiation to ingoing 
radiation is to be expected. However, the numerical miracle discussed in Sect. IV indicates 
that the backscattered ingoing radiation vanishes exactly at the future similarity horizon 
f = — 1, for reasons that we don't understand. In any case, an irregular pulse of outgoing 
radiation propagates outwards, with an infinite number of wavefronts piling up at the 
future similarity horizon. This is the nature of the pulse that a distant observer sees. 

This irregular pulse originates in the "echoing" behavior of the collapse in Region I. 
Choptuik found the "echoing" scale-factor 

30 

in his real choptuon, Eqs. (2). In the complex choptuon the echoing behavior is somewhat 
different, but the comparable scale for the solution to recur is measured by 2'n radians in 
the phase oscillation of so that 

gA ^ g27r/a; ^ 26.583086 ± 0.000005, 
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a value which clearly differs from Choptuik's. Therefore, this critical exponent A is not 
universal between the two solutions. 

There is an exact self-similar solution for a collapsing pulse of real scalar field in 
spherically symmetric general relativity [17,3,4], found independently by at least three 
groups. It would not be surprising to learn of further independent discoveries. This 
solution — which we shall call the MBONT solution, after its several discovers — exhibits 
critical behavior of the black hole mass at threshold [3,4]. However, the MBONT solution 
appears different from the choptuons in important respects. It does possess a past similarity 
horizon; however the scalar field ^ is not regular but irregular there — having a step in 
its derivative. A closely related issue is that the MBONT solution does not evolve from 
regular initial conditions for gravitational collapse. It remains to be seen whether or not 
the MBONT solution is an attractor. We shall return to this solution in a future paper. 

In this paper we have worked exactly at the threshold for black hole formation, and have 
not addressed the critical behavior of the black hole mass, Eq. (1). This critical behavior 
can be determined by first order perturbation theory around the complex choptuon, just 
as it can be addressed for the Evans and Coleman hot-gas choptuons [11,18]. Perturbation 
theory can also determine whether or not the solution is an attractor. We will return to 
these issues in a future paper. 

Region III, to the future of the singular point, is exactly fiat in our solution. Since the 
future similarity horizon is a Cauchy horizon, evolution cannot be unique; according to 
classical relativity, anything could come out of the singular point. However, it is remarkable 
that the solution admits an evolution into exactly fiat spacetime, and we have argued that 
this is the preferred evolution. 

Thus our solution provides one possible answer to the question: What is the endstate 
of gravitational collapse at the threshold for black hole production? The endstate is an 
outgoing, irregular pulse of radiation from the singular point in the wave zone , together 
with the most conservative of all possible endstates in the near zone: Flat spacetime. 

Further physical insight about this question must ultimately come from quantum grav- 
ity or string theory, since strong curvatures comparable to the Planck value occur near the 
singular point. 
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APPENDIX A 



In Table 1 we present all numerical solutions returned by our algorithm. As explained 
in Sect. IV, accuracy is lost near the past similarity horizon z = z^-, due to decay of 
the irregular part of (/> there, Eqs. (23). Physically we would like the solution to be 
at 2 = ^2- However, our numerical algorithm cannot distinguish between degree 
of smoothness and C°°, because Taylor expansions of the regular solution around 
were carried out only through second order. (If we had carried them to n*^ order, the 
algorithm could distinguish from C°°, but could not distinguish C"^"*^.) In turn, the 
degree of smoothness of the irregular solution depends on Eqs. (23), in such a way that 
our algorithm is reliable for /t<l, but loses reliability for k 1/3. The solutions marked 
"Reliable? No" in Table 1 all have k 1/3 and therefore we do not know whether is 
irregular and C"^, or regular and C°°, at 2; = 2:2 for them. The "Interpretation" depends on 
the behavior of the solution outside the past similarity horizon at 2; = 2:2. Solutions marked 
"Black Hole" encounter a numerical singularity that appears to be an apparent horizon. 
Such an apparent horizon probably means that the solution represents, not a choptuon, 
but a black hole that is growing by self-similar accretion of scalar field. Solutions marked 
"Choptuon" encounter no such apparent horizon, but do encounter a future similarity 
horizon. 



TABLE 1. Numerical solutions found. Solution is started at 2; = with values u) 
and P, in Eqs. (25). The past similarity horizon appears at 2 = 22- All solutions 
marked "Reliable? No" have k 1/3 there, and hence may have lost accuracy. 





P 


^2 


Reliable? 


InterpTetaiion 


0.0 


0.16286 


1.10668 


Yes 


Robertson- Walker solution 


0.0 




1.24793 


No 


Black Hole? 


0.15265 


0.23468 


1.35664 


No 


Black Hole? 


0.29825 


0.24502 


1.40583 


No 


Black Hole? 


0.43076 


0.25919 


1.47793 


No 


Black Hole? 


0.55509 


0.27670 


1.57717 


No 


Black Hole? 
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0.68307 


0.29923 


1.72316 


No 


Black Hole? 


0.81757 


0.32998 


1.97453 


No 


Black Hole? 


0.94825 


0.37269 


2.49240 


No 


Choptuon? 


1.04980 


0.44740 


4.56352 


No 


Choptuon? 


1.91544 


0.67217 


5.00354 


Yes 


Complex Choptuon; 



REFERENCES 

1. M.W. Choptuik, Phys. Rev. Lett. 70, 9-12 (1993). 

2. A.M. Abrahams & C.R. Evans, Phys. Rev. Lett. 70, 2980-2983 (1993). 

3. RR. Brady, Classical and Quantum Gravity 11, 1255-1260 (1994). 

4. Y. Oshiro, K. Nakamura & A. Tomimatsu, Progr. Theo. Phys. 91, 1265-1270 (1994). 

5. A. Strominger & L. Thorlacius, Phys. Rev. Lett. 72, 1584-1587 (1994). 

6. J. Traschen, Discrete Self-Similarity and Critical Point Behavior in Fluctuations 
about Extremal Black Holes, preprint ^r-qc/ 94030 16| , 1994. 



7. Garfinkle, D. Choptuik scaling in null coordinates, preprint |gr-qc/9412008| , 1994. 



8. S.W. Hawking & G.F.R. Ellis, The Large Scale Structure of Space- Time (Cambridge 
University Press, Cambridge, 1973). 

9. Penrose, R. 1979. In General Relativity, an Einstein Centennary Survey. S.W. Hawk- 
ing and W. Israel, Eds.:581-638. Cambridge University Press. Cambridge. 

10. See also D. Christodoulou, Comm. Math. Phys. 105, 337-361 (1986); 106, 587-621 
(1986); 109, 591-611, (1987); 109, 613-647, (1987). 

11. C.R. Evans, & J.S. Coleman, Phys. Rev. Lett. 72, 1782-1785 (1994). 

12. M. Choptuik, private communication (1994). 

13. D.M. Eardley, Comm. Math. Phys. 37, 287 (1974). 

14. B.K. Berger, J. Math. Phys. 17, 1268 (1976). 

15. D. Eardley, J. Isenberg, J. Marsden & V. Moncrief, nomothetic and Conformal Sym- 
metries of Solutions to Einstein's Equations, Comm. Math. Phys. 106, 137-158 
(1986). 

16. W.H. Press, B.P. Flannery, S.A. Teukolsky, and W.T. Vetterling, Numerical Recipes 
(Cambridge University Press, Cambridge, 1986). 

18 



17. Maithreyan, T., unpublished Ph.D. Thesis, Boston University, 1984. 

18. C. Evans, private communication (1994). 



19 



FIGURE CAPTIONS 



Figure 1. Coordinate systems in Minkowski spacetime. Arrows show the homothetic 
kiUing vector field ^ = tdt + rdr- Shown also are the past and future light cones of the 
origin; in the terminology of this paper, these light cones are past and future similarity 
horizons, where ^ becomes null. For coordinates z and v, see Sect. III. 

Figure 2. Behavior of the complex choptuon, the solution for oj = 1.9154446, over the 
domain < 2; < 10. The past similarity horizon is located at Z2 = 5.0035380, and oj was 
determined by demanding regularity there, (a) The complex functions q{z) and p{z) which 
represent the scalar field 0. (b) The metric functions b{z) and u{z). 

Figure 3. Behavior of the complex choptuon, the solution for uj = 1.9154446, over the 
domain —0.8 < < 10. The future similarity horizon is located at the value 2:4 of z where 
v{z) = —1. (a) The complex functions Q{w), P{w) which represent the scalar field (j). 
They oscillate infinitely many times in approaching the future similarity horizon, (b) The 
metric functions v{z) and u{z). Note that u{z^) = to numerical accuracy, showing that 
the solution matches onto flat spacetime at the future similarity horizon. 

Figure 4. Same as Fig. 3, except that the horizontal axis is now logarithmic, and shows 
ln(l + (a) The complex functions Q{w), P{w) which represent the scalar field (p. They 
oscillate infinitely many times in approaching the future similarity horizon, (b) The metric 
functions v{z) and u{z), plotted as ln(l + v) and In('u), both of which approach on the 
future similarity horizon. 

Figure 5. Interpretation of the complex choptuon. The similarity horizons lie at 
V = +1 (past) and v = —1 (future). Dotted lines show peaks and valleys of the scalar 
field {q for v < —1, Q for v > —1), to illustrate its oscillations. Region I is a collapsing 
sphere of gravitationally bound scalar field; in Regions Ila and lib this blends smoothly 
into an outgoing scalar wave. The outgoing scalar wave oscillates infinitely many times 
approaching the future similarity horizon at f = — 1. Region III is fiat. 
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